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Abstract 

Lattice-free sets (convex subsets of without interior integer points) and their 
applications for cutting-plane methods in mixed-integer optimization have been stud¬ 
ied in recent literature. Notably, the family of all integral lattice-free polyhedra 
which are not properly contained in another integral lattice-free polyhedron has 
been of particular interest. We call these polyhedra Z'^-maximal. 

It is known that, for fixed d, the family Z'^-maximal integral lattice-free polyhedra 
is finite up to unimodular equivalence. In view of possible applications in cutting- 
plane theory, one would like to have a classification of this family. However, this 
turns out to be a challenging task already for small dimensions. 

In contrast, the subfamily of all integral lattice-free polyhedra which are not 
properly contained in any other lattice-free set, which we call R'^-maximal lattice- 
free polyhedra, allow a rather simple geometric characterization. Hence, the question 
was raised for which dimensions the notions of Z’^-maximality and R.'^-maximality 
are equivalent. This was known to be the case for dimensions one and two. On the 
other hand, Nill and Ziegler (2011) showed that for dimension d > 4, there exist 
polyhedra which are Z'^-maximal but not R^-maximal. In this article, we consider 
the remaining case d = 3 and prove that for integral polyhedra the notions of R^- 
maximality and Z^-maximality are equivalent. As a consequence, the classification 
of all R^-maximal integral polyhedra by Averkov, Wagner and Weismantel (2011) 
contains all Z^-maximal integral polyhedra. 
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1 Introduction 

We call a d-dimensional convex subset C of lattice-free if the interior of C contains 

no points of For a subset X of we call a lattice-free convex set C X-maximal 
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if for every x £ X \ C, the set conv(0 U {x}) is not lattice-free. For Y C we say 
that P is a y-polyhedron if P = conv(Pny). A Z'^-polyhedron is also called an integral 
polyhedron. We denote the family of all P-polyhedra by ViY). In particular, P(R'^) is 
the family of all polyhedra in R"^ and is the family of all integral polyhedra in R'^. 

Recently, the family of Z'^-maximal integral lattice-free polyhedra has attracted at¬ 
tention of experts in algebraic geometry and optimization; see [TreOSllTrein] . [NZll] and 
|AWW11| . With a view towards applications in the mentioned research areas, having 
a better geometric understanding of such polyhedra is desirable. Currently, we have a 
rather good geometric description of R'^-maximality (which is a stronger property than 
Z'^-maximality) for general lattice-free sets; see [Lov 89] and [A vel3| . In contrast, so far 
no simple description of Z'^-maximality is available, even in the case of integral polyhedra. 
This has motivated the following question, asked in [NZllj and [AWWlI] : is every Z'^- 
maximal integral lattice-free polyhedron also R'^-maximal? This can also be formulated 
as follows: is it true for every integral lattice-free polyhedron P, that if conv(PU {x}) is 
lattice-free for some x € R'^ \ P, then conv(P U {x}) is lattice-free for some x G Z'^ \ P? 
The answer to this question is trivially ‘yes’ for dimensions d = 1,2, while it was shown 
in [NZllj that the answer is ‘no’ for each d > 4. The purpose of this paper is to settle 
the case d = 3, which remained open; see [NZlll Question 1.5] and [Wagll Section 7.4]. 


Theorem 1. Every -maximal integral lattiee-free polyhedron is also -maximal. 


The interest in the family of lattice-free integral polyhedra, and especially those 
that are Z^^-maximal, was raised by potential applications in mixed-integer optimization, 
more precisely in cutting-plane theory. It is known that lattice-free sets can be used in 
mixed-integer programming for describing various families of cutting planes. Consider 
the set P n X R”) of feasible solutions of an arbitrary mixed-integer linear program, 
where P C is a rational polyhedron, d G N is the number of integer variables and 

n G N U {0} is the number of continuous variables. Given a lattice-free set L C R'’^, we 
call a closed halfspace H of an L-eut for the polyhedron P if P \ int(P) x R"" C H. 
Furthermore, given a family L of lattice-free subsets of R'’*, we say that P is an L-eut if 
H is an P-cut for some L £ L. 

In the theory of general mixed-integer linear programs, one is interested in explicit 
descriptions of families L being possibly simple and small on the one hand, and generating 
strong cuts on the other hand. The strength of £-cuts can be represented in terms of 
the sequence of P-closures, defined as follows. For a rational polyhedron P C RP*^, we 
define the L-closure of P by cc{P) = (Q C P : P is an P-cut}. For a nonnegative 
integer k, the k-th L-elosure c^{P) of P is defined recursively by c^{P) = cjr{c^^(P)) 
for A: G N and c^{P) = P. Clearly, c^{P) contains the mixed-integer hull Pmi := 
conv(P n (Z'’* X R”)) of P as a subset. 

The strength of P-cuts can be formalized by describing how well Pmi is approximated 
by the sequence (c£(P))^g. The strongest form of approximation is expressed as finite 
convergence to the mixed-integer hull. In the case of finite convergence, c^{P) coincides 
with Pmi if k is large enough (where the choice of k depends on P in general). If the 
latter holds for every rational polyhedron P C RP**, we say that P has finite convergence 
property. Note that, a priori, the finite convergence property depends not only on P but 
also on the number of continuous variables re. 
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As an example, the family of split sets that corresponds to the classical split cuts 
has finite convergence property for n = 0 (this follows from the fact that the well-known 
Gomory cuts are a subclass of split cuts and one has finite convergence for Gomory cuts; 
see e.g. [Sch86t §23]) but does not have finite convergence property for every re > 0 (see 
|GKS9nj l. Thus, more complicated lattice-free sets have to be used for achieving finite 
convergence in the mixed-integer setting. Del Pia and Weismantel showed in |DW12j that 
the family of all integral lattice-free polyhedra has finite convergence property for every re. 
Since, for a general d, the family of all integral lattice-free polyhedra is complicated and 
hard to describe explicitly, having a smaller family of lattice-free polyhedra with finite 
convergence property for every re is desirable. It turns out that every integral lattice-free 
polyhedron is a subset of a Z'^-maximal integral lattice-free polyhedron. This is stated 
without proof in |DW12j and can be proven using results in |NZ11] or arguments from 
[ AWW1 f| . Gonsequently, the family of all Z'^-maximal integral lattice-free polyhedra 
has finite convergence property for every re. Furthermore, no proper subfamily of the 
latter family has finite convergence when re > 0, as follows from results of |Dell2j . In 
view of the above comments, it would be interesting to have an explicit description of 
all Z'^-maximal integral lattice-free polyhedra. 

In |AWWlf| and [NZllj it was shown that, there exist finitely many polyhedra 
Li,..., Ln with N depending on d such that every Z'^-maximal integral lattice-free 
polyhedron L is an image (f){Li) of some Lj, zG {!,...,A^}, under an affine transforma¬ 
tion (j) satisfying (/'(Z'^) = Z'^. Affine transformations 4> of satisfying 0(Z'^) = Z'^ are 
called unimodular, while subsets A, B of satisfying (()(A) = B for some unimodular 
transformation (p are said to be unimodularly equivalent (we write A B for sets A, B 
which are unimodularly equivalent). Thus, the family of Z'^-maximal integral lattice- 
free polyhedra is finite up to unimodular equivalence. The problem of enumerating this 
(essentially finite) family for small dimensions was raised in [AWWlf] . For dimensions 
d = 1 and d = 2 enumeration can be carried out easily. Already, for dimension d = 3, 
this is a challenging problem which was only partially addressed in [AWWlI] . Using 
properties of M'^-maximal lattice-free sets, the authors of |AWWlI] were able to enumer¬ 
ate the family of all M^-maximal integral lattice-free polyhedra. While, a priori, this is 
only a subfamily of the family of Z^-maximal integral lattice-free polyhedra, our Theo¬ 
rem [H shows that the classification given in [AWWllj also provides an enumeration of 
Z^-maximal integral lattice-free polyhedra. 

Corollary 2. Up to unimodular equivalence, there are exactly twelve bounded Z^- 
maximal integral lattice-free polyhedra (see Figures\^ andWi) o,nd two unbounded ones, 
namely [0,1] x and conv((0,0), (2,0), (0,2)) x M. 

We remark that the main property of M'^-maximal lattice-free polyhedra which was 
exploited in the classification provided in [AWWlf] is the following: a lattice-free poly¬ 
hedron P is M'^-maximal if and only if each of its facets contains an integral point in its 
relative interior; see |Lov89] . We call a facet with the latter property blocked and we will 
use this characterization of R'^-maximality repeatedly throughout this article. 
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Figure 1: The Z^-maximal integral lattice-free polytopes with lattice width two. For 
further reference, the polytopes are labeled by a pair of indices where i is the 

number of facets and j the lattice diameter (defined at the end of the introduction). 



Figure 2: The Z^-maximal integral lattice-free polytopes with lattice width three. 

Proof strategy 

In the proof of Theorem[Tl we use a classification of all Z^-maximal polytopes in P(^Z'^). 
This is provided in Section [2l Every such polytope is contained in an M^-maximal lattice- 
free convex set L in the plane and its vertices then have to be contained in L n We 
give a slightly extended version of the well-known classification of M^-maximal lattice-free 
convex sets L which allows us to enumerate all Z^-maximal lattice-free ^Z^-polyhedra. 
We then turn to integral Z^-maximal lattice-free polyhedra in dimension three. We 
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can restrict ourselves to polytopes here since every unbounded Z^-maximal lattice-free 
polyhedron in P(Z^) is unimodularly equivalent to [0,1] x or conv(o, 2ei, 2e2) x R 
in view of [AWWlIl Proposition 3.1], The R^-maximality of those two sets is easy to 
see. Furthermore, one can easily observe that a lattice-free polytope in P(Z^) with 
lattice width one (see page ED cannot be Z^-maximal, since any such polytope is prop¬ 
erly contained in an unbounded lattice-free integral polyhedron formed by two adjacent 
parallel lattice hyperplanes of Z^. Therefore, it suffices to consider only the polytopes 
with lattice width at least two, since the lattice width of an integral polytope is inte¬ 
gral. We distinguish two cases: polytopes with lattice width two and polytopes with 
lattice width at least three. In the case of lattice width two, we first analyze Z'^-maximal 
lattice-free integral polytopes with this lattice width in general dimension. Assuming a 
Z'^-maximal lattice-free polytope P to be contained in R'^"^ x [—1,1], we deduce prop¬ 
erties of Po = {x ^ R*^"^ : (x,0) G P} and on the faces of P contained in R^ x {—1} 
and R^ x {!}, respectively. In particular, we obtain that Pq is a Z'^“^-maximal ^Z'^“^- 
polytope. Narrowing the focus to the case of dimension three, we then make use of the 
fact that Pq has to be one of the four Z^-maximal lattice-free ^Z^-polytopes which we 
enumerated in Section [2l Based on that, we classify all Z^-maximal polytopes in P(Z^) 
which have lattice width two. It turns out that every such polytope is also R^-maximal. 

In Section 01 we complete the proof of Theorem [T] by a computer search. To this 
end, we develop an algorithm which finds all Z^-maximal polytopes in P(Z^) with lattice 
width at least three and verifies that all of them are also R^-maximal. For enumerating 
all such polytopes, we establish bounds on the volume, the first successive minimum and 
the lattice diameter of Z^-maximal integral lattice-free polytopes with lattice width at 
least three, making use of classical results from the geometry of numbers. Finally, an 
implementation of the search algorithm is provided in the appendix. 

Basic notation and terminology 

For background information on the theory of mixed-integer optimization, convex sets 
and polyhedra, and the geometry of numbers, see for example |Sch86| . [Bar 02], [R PC 97], 
|Sc.h93] and [CTT^ . 

Throughout the text, d G N is the dimension of the ambient space R'^, which is 
equipped with the standard Euclidian scalar product, denoted by (•, •). By o we denote 
the zero vector and e, denotes the i-th unit vector, where in both cases the dimension 
is given by the context. For two points x, y G R'^, we denote the closed line segment 
connecting those points by [x,y]. Accordingly, one has the open line segment {x,y) := 
{Ax -|- (1 — A)y : 0 < A < 1}. For a subset X of R'^, we write Y C A if E is a subset 
of X and Y ^ A if, furthermore, Y ^ X holds. The cardinality of A is denoted by 
|A|. By —A we denote the set {—x : x G A}. For two subsets X,Y of R'^, we write 
A-|-E :={x-|-y : xGA, yGE} for the Minkowski sum of A and Y, while A — E is 
short for A -|- (—E). For t G R, we write tX := {tx : x G A}. For a set Z C Z, we 
denote by gcd(Z) the greatest common divisor of the set Z. We call a vector z G Z‘^\{o} 
primitive if the greatest common divisor of the components of z is one. 

For a set A C R"^, we use conv(A), aff(A) and hn(A) to denote the convex hull, 
affine hull and linear hull of A, respectively. By bd(A), relbd(A), rehnt(A) and int(A) 
we denote the boundary, the relative boundary, the relative interior and the interior 
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of a convex set K, respectively. By vol(ii'), we denote the volume of K and dim(ii') 
denotes the dimension of K, which is defined as the dimension of the affine hull of K. A 
convex set K such that K = —K is said to be o-symmetric. A compact convex set with 
non-empty interior is called convex body. For a convex body K C the convex body 
K + {—K) is called the differenee body of K. For a convex set iF C we say that K' 
is a homothetic copy of K K' = XK -|- v for some A G M and v G We call K' a 
positive (non-negative) homothetic copy of iF if A is a positive (non-negative) number. 
For a convex set iF C and a vector u G \ {o}, we denote by h{K,u) the support 
function h{K,u) := sup{(u, x) : x G K}. 

Let bi,...,bd linearly independent elements of Then we call the group A = 
{ 2:161 Zdbd ; 2 : 1 , ..., 2 ;^ G Z} a lattice of rank d. The vectors 61 ,..., 6 ^ are said to 

form a lattice basis. The determinant of A is defined as the determinant of the matrix 
having columns 61 ,..., 6 ^. We say that a subgroup iL of a lattice A is a lattiee hyperplane 
if H spans a {d — l)-dimensional affine subspace of Two parallel lattice hyperplanes 
are said to be adjaeent if there are no points of A properly between them. 

For a subset X of and a vector u G by 

w(JF, u) = sup {u , x) — inf {u , y) 
xex 

we denote the width funetion of X. If u is a primitive integral vector, then w(X, u) is 
called the (lattiee) width of X in the direction of u. Note that if X is compact (say, a 
convex body), the infimum and supremum, respectively, are attained. The lattiee width 
of a set X C is defined as the infimum of the width function of X over all non-zero 
integral vectors, i.e. 

lw(A) := inf w{X,z). 

Again, for a convex body the infimum is attained. We say that the lattice width of a 
convex body K is attained for z' gU^\ { 0 } if Iw(iF) = w(iF, z). 

A polyhedron P C is the intersection of finitely many closed halfspaces. A 
bounded polyhedron is called a polytope and a two-dimensional polytope is called polygon. 
A polyhedron P is called rational if there exists s G N such that P G that is, 

P = conv(P n ^Z'^). 

For a polyhedron P and a vector u G M'^\{o}, if h{P, u) is hnite, we dehne F(P, u) := 
{x G P : {x , u) = h{P,u)} and call F{P,u) a face of P. If F := F{P,u) has dimension 
dim(P) — 1, we say that F is a facet of P and say that u is an outer normal (facet) 
vector of F. If F consists only of a single point, then this point is called a vertex of P. 
The set of vertices of P is denoted by vert(P). For a d-dimensional rational polyhedron 
P, we define U (P) to be the set of all primitive vectors which are outer normal vectors 
of facets of P. Note that if P is rational, every facet of P can be given as F{P,u) 
for some u G U{P). Let x G bd(P). Then the normal cone of P at x is dehned as 
N{P,x) := { 0 } U {u G \ { 0 } : x G P(P, n)}. We call a vertex u of P unimodular 
in this article if v has the property that there are exactly d facets of P containing v 
and these facets can be given as P(P, ui),... ,F{P, Ud), where ui,... ,Ud are elements of 
U (P) forming a basis of the lattice Z^. 

A d-dimensional polytope with d-|-1 vertices is called simplex of dimension d. Given 
a d-dimensional simplex S' C with vertices vi,... ,Vd+i 3 -nd a point x G R*^, the 
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barycentric coordinates of x with respect to S are uniquely determined real numbers 
Pi,..., Pd+i satisfying x = Yjt=i Pi = 1- case, for i G {1 ,..., d+l}, 

we say that Pi is the barycentric coordinate associated with Vi. One has x G int(5') if 
and only if /3i > 0 for all i G {1,..., d + 1}. 


2 Half-integral lattice-free polyhedra in dimension two 

As explained in the outline of the proof, in this section we show that for two-dimensional 
half-integral lattice-free polygons, Z^-maximality implies M^-maximality. 

We make use of the following characterization of M^-maximal lattice-free convex sets 
in Although the characterization does not seem to be available in the literature in 
precisely this form, it is largely known from |Hur 90 j ■ [DWlOl Proposition 1 ] and |AW 12 j . 

Theorem 3 (Classification of maximal lattice-free sets in dimension two). Let L he an 
^-maximal lattice-free convex set in Then there exists a unimodular transformation 
(f such that ip{L) G £i U ... UTs, where Ci,... ,C^ are the following families (see Fig. 

(a) Cl consists of exactly one polyhedron, namely the unbounded set [0,1] x R. 

(b) £2 is the set of all ‘axis-aligned’ triangles with vertices ui,U2,U3 such that vi coin¬ 
cides with the vertex o of the square [0,1]^, the vertices V2,V3 satisfy 

V2 G (l,+oo) X { 0 }, Vs G { 0 } X (l,-|-oo) 

and for the vertices ei, 62, ei -|- 62 of the square [0,1]^ one has 

eie{vi,V 2 ), e2G(ui,U3), ei62 G (u2, U3). 

(c) £3 is the set of triangles with vertices vi,V2,vs satisfying 

vi G (-00, 0 ) X { 0 }, V2 G ( 1 , -1-00) X { 0 }, Vs G ( 0 , 1 ) x ( 1 , -|-oo) 
and such that for the vertices o, ei, 62, ei -|- 62 of the square [0,1]^ one has 
o,ei G (ui,U2), e2G(ui,i;3), ei-|-62 G (u2, U3). 

(d) £4 is the set of triangles with vertices vi,V2,vs satisfying 

vi G [{xi,X 2 ) G R^ : X 2 > 1 , 0 < xi -I- X 2 < 1} , 

V2 G {ixi,X2) G R^ : 0 < Xi < 1 , Xi -t- X2 < 0} , 

Vs € {ixi,X2) G R^ : xi > 1 , 0 < X2 < 1} 

and such that for the vertices o, ei, 62 of the triangle conv(o, ei, 62) one has 

o G {vi,V2), ei G iv2,vs), 62 G (u3,ui). 
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(e) £5 is the set of quadrilaterals with vertices fi, ^^ 3 , ^^4 satisfying 


vi G (-(X),0) X (0,1), V 2 G (0,1) X (- 00 , 0 ), 

V 3 G (1, + 00 ) X (0,1), V 4 G (0,1) X (1, + 00 ), 


and such that for the vertices o, ei, 62 , ei + 62 of the square [ 0 , 1 ]^ one has 
oe{vi,V2), eie{v2,V3), ei + 62 G (i; 3,7;4), e2G(t>4,fi). 

£2 


0000 



Vl V2 




Figure 3: Illustration of the classes £ 2 , • • •, £5 of Theorem [3l The area shaded in darker 
grey is [0,1]^ for £ 2 , £ 3 , £5 and conv(o, ei, 62 ) for £ 4 . The areas shaded in light grey 
are the regions containing the vertices of </?(£)• 


Proof. If L is unbounded, then it is well known from [Lov891 §3], it is unimodularly 
equivalent to [0,1] x M G £ 1 . If £ is bounded, from [Lov89l §3] we have that L is 
either a triangle or a quadrilateral and again from [Lov891 §3] we have that each edge 
of L is blocked. It is well-known (see, e.g., |Hur90j . |DW10l Proposition Ij) that if L 
is a quadrilateral, then the four integral points blocking its edges can be unimodularly 
mapped onto {0,1}^. It is then easy to see that this mapping transforms L into an 
element of £ 5 . This leaves the case that L is a triangle. Assume first that L has an 
integral vertex v. Then we can choose an integral point in the relative interior of each of 
the three edges of L such that those three points together with v can be unimodularly 














mapped onto { 0 , 1 }^ (see again |DW 10 [ Proposition 1 ]) where v is mapped onto o. Thus, 
(/?(L) G £2- We now switch to the case that L does not have integral vertices. If 
L contains at least four integral points, then the relative interior of some edge of L 
contains at least two integral points ui,U2 such that (tti,tt2) H = 0 . Together with 
one point from the relative interior of each of the two remaining edges, they form an 
integral quadrilateral not containing any other integral points. Unimodularly mapping 
this quadrilateral onto {0,1}^ such that ui and U2 are mapped onto o and ei, respectively, 
yields an element of £3. 

This leaves the case that L has precisely three integral points. There exists a unimod- 
ular transformation cp such that (p{L) n Z^ = {0,61,62} and (p(L) has vertices vi,V2,V3 
with o G (^1,^2), 61 G {v2,V3) and 62 G (n3,ni) (see, e.g., |Hur 90 ] l. We prove p{L) G £4, 
for which it remains to show that fi, ^2, ^3 are situated as claimed in |(d)[ We can express 
0,61,62 with respect to vi,V2,V3 as follows: there exist 0 < ai,0:2,03 < 1 such that 

o = aivi + (1 - ai)v2, 

61 = a2V2 + (1 - a2)v3, (1) 

62 = asvs + (1 - a3)vi. 


We now want to express vi,V2,V3 with respect to 0 , 61 , 62 . We follow the argumentation 
in [AW 12 ( Proof of Lemma 4 . 2 ] and [Hnr 9 n( ( 7 )] and translate ([I]) into the following 
matrix representation: 




CKi 0 1 — 03^ 

1 — cx\ 0.2 0 

0 1 — 0:2 0:3 


The determinant of the matrix involving ai, a 2 , 0:3 is D := aia 2 CK 3 + (1 — a:i)(l — q: 2)(1 — 
03 ) > 0. Hence, this matrix is invertible and we can directly compute the components 
of vi,V 2 -,V 3 via Cramer’s rule: 

t’l = - (1 - ai)a3, (1 - ai)(l - 02 )), 

6’2 = ^ (aia 3 , -ai{l - 02 )) , ( 2 ) 

t’s = ^ ((1 - ai)(l - as), aia 2 ) • 

In [AW12t Lemma 5.1] it was shown that either 


Ui + aj <1 for all 1 < i < j < 3 (3) 

or ai + aj > 1 for all 1 < i < j < 3. In the case that dS]) holds, it is straightforward to 
verify that the vertices vi,V 2 ,V 3 are situated as claimed. If a* aj > 1 for every i,j G 
{1,2,3} holds instead of ([3]), one can instead prove in a straightforward way that this 
yields the following regions for the vertices vi,V2,V3: 

vi G {(xi,X2) G : xi + X 2 < 0 , 0 < X 2 < 1} , 

V 2 £ {{xi,X 2 ) G : xi > 1 , 0 < xi + X 2 < 1} 

V3 G {(xi,X 2 ) G : 0 < xi < 1 , X 2 > 1} . 
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Applying the unimodular transformation {xi,X 2 ) {x 2 ,xi), and appropriately re¬ 

indexing the vertices yields a triangle in £ 4 . This is again straightforward to verify. □ 

Based on the previous theorem, we now enumerate all Z^-maximal lattice-free half¬ 
integral polygons up to unimodular equivalence. The polygons Q 2 ,Q 3 ,Q 4 ,Q 5 of the 
following theorem are depicted in Figure 01 



Figure 4: All Z^-maximal lattice-free ^Z^-polytopes. 


Theorem 4 (Z^-maximal lattice-free polyhedra in P(^Z^)). Let Q G P(^Z^). Then Q 
is 1?-maximal lattice-free if and only if Q Qi for some i G {1,..., 5}, where 

Qi is the unbounded polyhedron [0,1] x M, 

Q 2 is the triangle with vertices (0,0), (2,0), (0, 2), 

Qs is the triangle with vertices (|±|, 0 ),(^,|), 

Q 4 is the triangle with vertices ± 1 , 0 ), (|, 2 ), 

Qs is the quadrilateral with vertices (— 5 , 5 ), (^, |), (|, 5 ), 

In particular, Q is Tp'-maximal lattice-free if and only if Q is M.'^-maximal lattice-free. 

Proof. Let Q G be Z^-maximal lattice-free. We show that then, Q '24 for 

some i G {1,... ,5}. It is known, see for example [AW121 Proposition 3.1], that every 
lattice-free set in is a subset of an R^-maximal lattice-free set. Thus, there exists 
an R^-maximal lattice-free set L such that Q C. L. By applying a suitable unimodular 
transformation, we can assume L G Ci with i G {1,. .., 5}, where £ 1 , ..., £5 are defined 
as in Theorem [3l We observe that one has 

conv(L n Z^) C Q (4) 

and 

Q C conv (^L n , (5) 

where (jl]) holds since Q is Z^-maximal and ([5]) holds since Q G P(^Z^) and Q Q L. 
Case 1: £ G £ 4 . Then £ = [0, 1] x R and hence, (0]) together with (0]) yields Q = L. 
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Case 3 


Case 2 




Case 5 


Figure 5: Illustration of the proof of Theorem HI Q is contained in the shaded area, the 
white dots are the elements of and the black dots are the elements of ^1? involved 
in the proof. 


Case 2: L £ € 2 - In view of (jl]), we have [0,1]^ C Q. Prom the definition of C 2 , we 
have that every point of L n (|Z^ \ [0,1]^) has one of the following forms: 

(3/2,0) +A: (1/2,0), (0, 3/2) + A: (0,1/2) , 

(3/2,1/2) + k (1/2,0), (1/2,3/2) + k (0,1/2) , 

where k is some non-negative integer; see also Figure [5l Let us first assume (|, |) and 
(|, |) are both in Q. Since (1,1) is contained in the relative interior of the line segment 
between these two points, [(|) 5 )) (^j |)] is contained in an edge of Q. Consequently, 
Q C conv ((0,0), (2,0), (0, 2)) = Q 2 and since Q is Z^-maximal, this implies the equality 
Q = Q 2 - Let us now assume that (^,1) or (|,|) is not in Q. By applying a linear 
unimodular transformation that swaps ei and 62 if necessary, we can assume (^, |) 0 Q- 
This implies that no point of the form (^, |) + A: (O, |) is in Q. Likewise, since (1, 1) £ Q 
and (|, |) G [(1,1), (0, 2)], we have (0, 2) ^ Q and consequently, (O, |) + A: (O, |) ^ Q for 
every A; > 1. Hence, either (0,3/2) £ Q or Q C [0,+oo) x [0,1]. If (0,3/2) £ Q, we have 
Q C conv ((0,0), (0,3/2), (3,0)) ~ Qs and the Z^-maximality of Q implies the equality 
Q ~ Q 3 . Otherwise, we then have Q C [0, +00) x [0,1] ^ M x [0,1], a contradiction to 
the Z^-maximality of Q. 

Case 3: L £ C^. In view of 0 , we have [0,1]^ C Q. From the definition of T 3 , we 
have that every point in L n (^Z^ \ [0,1]^) has one of the following forms: 

(3/2,0) +A: (1/2,0), (-1/2, 0) + A: (-1/2, 0) , (1/2,3/2) + A; (0,1/2) , 

(3/2,1/2) + k (1/2, 0) , (-1/2,1/2) + k (-1/2,0) , 
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where k is some non-negative integer; see also Figure O 

Analogously to Case 2 , we can now restrict the values of k for each of the five forms: 
if a point of any of the forms is in L for which the value of k is too large, then, because 
(0,1) and (1,1) are in the boundary of L, we have that L n ^7? is contained in M x [0,1] 
or [ 0 , 1 ] X M. In those cases, in view of ([5]), Q is contained in a lattice-free unbounded 
polyhedron in V{7?). Since Q is bounded, this inclusion is strict, a contradiction to the 
Z^-maximality of Q. This leaves a finite number of points in ^7? which can be in Q and 
from dS]), we obtain that the vertices of Q which are not in {0,1}^ are among the points 


/ X / lx/ 1 lx /I X /I 3, ,3 , ,3 1, , , 

( ~ 1)0)) ( “ QjO)) ( “ 2’ 2^’ (o’ 9)’ (o’®)’ (o’ o)’ (^’®)‘ 


- 2 ’ 2 ' 


' 2 ' 2 ' 


It can be checked in a straightforward way that the only Z^-maximal lattice-free polygons 
in with vertices in this set are Q 3 and Q 4 . 

Case 4- L G £ 4 . Here, we deduce A := conv ((0, 0), (0,1), (1, 0)) C Q from (jlj). 
From the definition of £ 4 , we have that every point of L n {\7? \ A) has one of the 
following forms: 


(1,1/2) +A: (1/2,0), (1/2,-1/2)+A: (0,-1/2), (-1/2,1) + A: (-1/2,1/2) , 


where k is some non-negative integer; see also Figure [5j As in the previous cases, we 
exploit the fact that (0,0), (0,1), (1,0) he in the boundary of L to obtain bounds on 
k for each of the three forms: if k is too large, we again deduce that Q is contained 
in an unbounded set unimodularly equivalent to M x [0,1]. Since Q is bounded, this 
contradicts the Z^-maximality oi Q. In fact, we have k <2 for each of the three forms. 
In view of ([5]), this yields a finite set of points of ^Z^ which can be vertices of Q. The 
straightforward verification that no lattice-free half-integral polygon with vertices in this 
set is Z^-maximal is left to the reader. 

Case 5: L £ C 5 . In view of 0 , we have [0,1]^ C L. Furthermore, it is easy to see 
that every point of {L n ^Z^) \ {0,1}^ has one of the following forms: 


1 1 

2 ’ 2 


± A: 


1 


,0 


1 1\ /I 

—, — I i A: ( 0, — 

2’2/ V 2 


where k is an integer and k > 2 ; see also Figure [5j One can observe that, since the 
points of {0,1}^ are in the boundary of Q, if there exists a point of any of the four 
forms with k > 2 we have that Q C £ n ^Z^ is contained in M x [0,1] or [0,1] x M. 
Since Q is bounded, this contradicts the Z^-maximality of Q. This yields Q ^ Q 5 . It 
is straightforward to check that every ^Z^-maximal integral polygon properly contained 
in Qs is not Z^-maximal, which implies Q = Qs- 

We have shown that every Z^-maximal polyhedron in P(^Z^) is unimodularly equiv¬ 
alent to Qi for some i G {1,..., 5}. It is easy to see that for every i G {1,..., 5}, the 
polyhedron Qi is also R^-maximal as all of its edges are blocked. This yields the equiv¬ 
alence of Z^-maximality and M^-maximality for lattice-free polyhedra in P(iZ^). □ 


Remark 5. One can ask the question whether Z'^-maximality and M'^-maximality are 
equivalent for lattice-free polyhedra in V{^7^) for given s G N. The previous theorem 
shows the equivalence for d = 2 and s = 2. This answers a question from [Wag 11 
Section 7.4]. Furthermore, since R(Z^) C R(iZ^), this equivalence also holds for d = 2 
and s = 1 . 
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3 Z^-maximal lattice-free integral polytopes with lattice 
width two 

We will now prove that every Z^-maximal lattice-free polyhedron in P(Z^) is also 
maximal. As explained in the proof strategy, we only need to prove this for polytopes, 
because for unbounded polyhedra in P(Z^) the equivalence of Z^-maximality and 
maximality is straightforward to see in view of [AWWill Proposition 3.1]. 

In the beginning of the previous section, we remarked that Z^-maximal lattice-free 
polytopes in P(^Z^) appear naturally when considering Z^-maximal lattice-free poly¬ 
topes in P(Z^). More precisely, they appear as hyperplane sections of Z^-maximal 
lattice-free polytopes in P(Z^), as we will explain in detail in this section. However, 
where possible, we do not restrict ourselves to dimension three but instead present the 
argumentation concerning lattice-free integral polytopes with lattice width two in ar¬ 
bitrary dimension d > 3. Our general setting is the following: let P G be a 

Z'^-maximal lattice-free polytope with lattice width two. Then there exists a primitive 
vector u G Z'^ \ {o} such that w(P, u) = 2. Changing the coordinates using an appro¬ 
priate linear unimodular transformation, we can always assume u to be and using 
some translation by an integral vector if necessary, we can furthermore without loss of 
generality assume P C x [—1,1]. 

Properties of faces and cross-sections 

If P G P(Z'^) is contained in X [—1,1], then the integral points of P are contained in 
]^d-i ^ |_2^o, 1} and in particular, the latter set contains all vertices of P. Note that if 
moreover P has lattice width two, then x {0} meets the interior of P and since P is 
lattice-free, this allows us to derive a characterization of Pr\W^~^ x {0}. We then provide 
relations between P n x {0} and F{P,ed),F{P,—ed)- For the characterization of 
P n X {0}, we need the following lemma. 

Lemma 6 . Let H be an affine subspace ofW^. Let A,B be convex subsets ofW^. Then 
the following statements hold: 

(a) If A Q H, then H n conv(A U B) = conv [Au {H Ci B)). 

(b) If A H, B % H and (1 — A) A -|- XB C H for some 0 < A < 1, then 

H n conv(A U P) = (1 — A)A -|- XB. 

Proof, (a): Let A F H. It suffices to verify the inclusion Pnconv(AuP) C conv (Au(Pn 
P)); the reverse inclusion is trivial, as both A and HCiB are subsets of H nconv(AuP). 
Consider an arbitrary x £ H Ci conv(A UP). That is, x is a point of H that can be 
written as a convex combination of some a £ A and b £ B. If x = a, then in particular 
X G A and we are done. If not, then b lies on the line passing through x and a. Since 
x,a £ H, this line lies in H and therefore, we also have b £ H. 

(b): Let A H and B H. It suffices to verify the inclusion H n conv(A U P) C 
(I — A)A -|- XB, as the reverse inclusion is trivial. Obviously, one has conv(A U P) = 
i]((f ~ + h^))- Now we fix an affine function / : —)• M such that H = 


13 




|x G : f{x) = O}. Since H / the function / is not identically equal to zero. 
Applying / to (1 — A)A + XB C H, we arrive at (1 — X)f{A) + Xf{B) C f{H), where 
by the choice of /, one has f{H) = {0}. In view of 0 < A < 1, the latter shows that 
f{A) and f{B) must be singletons. We interpret the singletons f{A), f{B), f{H) as real 
values, where f{H) = 0. Let f{A) = a and fiB) = /?. Since A and B are not subsets 
of H, we have a,/3 / 0 and in view of (1 — A)a + A/3 = 0, we have a ^ 13. Hence, for 
// G [0,1], one has (1 — fi)a + /r/3 = 0 if and only if ^ = A and, thus, we have shown that 
(1 — fi)A + fiB is disjoint with H unless /x = A. Consequently, we obtain (b). □ 

Theorem 7. Let P G V{h^) he a TJ^-maximal lattice-free polytope such that P C x 

[—1,1] and w{P,ed) =2. For / G {—1,0,1}, let Pi := {x ^ : (x,/)gF}. Then 

the following statements hold: 

(a) Pq belongs to and is ^-maximal lattice-free. 

(b) P-i and Pi belong to and have the following properties relative to Pq: 

(bl) Pi and F_i are integral polytopes satisfying 

Pi + F_i C 2Po. (6) 

(b2) The pair (Fi,F_i) satisfies the following maximality condition: for polytopes 
R-i,Ri G F(Z‘^“^), the inclusions Pi ^ Ri, P-i C R_i and Ri R-i C 2Fo 
imply Ri = Pi and R-i = P-i. 

(b3) Pi and F_i satisfy 

2conv( vert(Fo) \ Z'^“^) C Fi + P_i. (7) 

Proof. Throughout the proof, we will use the following notation: H denotes the hyper¬ 
plane X {0} and for i G {—1,0,1}, we write Pi := Pi x {/}. We will repeatedly 

make use of the fact that P C R'^"^ x [—1,1] and hence, F n Z*^ C Z'^"^ x {—1, 0,1}. In 
particular, taking into account P G F(Z'^), we have 

vert(F) C F_i U Fq U Fi. (8) 

The assumption w(F, e^) = 2 implies Fi / 0 and F_i / 0. 

Consider 


A := conv(i7 n vert(F)) and B := conv {P-i U Pi). 

From ([ 8 ]), it is clear that conv(A L) B) = P. Hence, in view of Lemmata), we get 

Pq = H D conv(A D B) = conv (A U (H n B )), (9) 

where the second equality is a consequence of Lemmata). 

Assertion (a): We hrst prove that Fq belongs to F(^Z'^), which is equivalent to 
vert(Fo) C x (Oj. In ([9]), the set A is an integral polytope, possibly empty. Thus, 

it suffices to verify that the vertices of HCiB belong to ^Z^“^ x {0}. In fact, Lemma[ 6 l[b) 
yields 

H n B = H D conv {P-i U Pi) = 
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Thus, since P_i and Pi are integral polytopes, P n P is half-integral and so is Pq. 

Next, we show that Pq is lattice-free. By [Roc971 Corollary 6 . 5 . 1 ] we have that 
relint(P D H) = int(P) Pi H and thus 

relint(Po) C = relint(P P P) P Z'^ = int(P) P P P Z"^ 

C int(P) P Z'^ = 0. 


Hence, relint(Po) does not contain points of Z*^"^ x {0} and thus, Pq is lattice-free. 

To complete the proof of |(a)[ it remains to show that Pq is Z‘^“^-maximal. Let 
p G X {0} \ Pq. As p ^ P and P is Z^-maximal lattice-free, there exists some 

z G in int (conv(P U {p})). Since P is lattice-free, we have z = Ap -|- (1 — X)x 
for some 0 < A < 1 and x G P. Let Xd,Pd,Zd denote the d-th component of x,p,z, 
respectively. Because Xd G [—1,1] and pd = 0, we obtain —1 < < 1 and hence Zd = 0. 

This implies 

z G relint (conv(P U {p})) P P = relint (conv(P U {p}) P P) 

= relint (conv(Po U {p})), 

where the first equality follows from |Roc971 Corollary 6.5.1] and the second equality 
follows from Lemma [6l] a). Hence, Pq is Z'^“^-maximal lattice-free. 

Assertion (b): Obviously, P_i,Pi are integral polytopes because P_i,Pi are faces 
of the integral polytope P. 

(bl): Clearly, Ap_i-|-iPi is a subset of both P and P. This yields ^P_i-|-^Pi C Pq 
and hence assertion (bl) follows. 

(b2): Let Pi,P_i C be integral polytopes such that Pi -|- P_i C 2Pq and 

Pi P Pi, P-i P P_i. We write Ri = Ri x {z} for i G {—1,1}. Consider the integral 
poly tope 

P := conv (P_i U Pi U Pq) , 

which is contained in R'^ 1 X [-1, 1] and contains P. In view of Lemma [61(a), we have 
P P P = conv ^Pq U conv ((Pi U P_i) P P)^ • 

Applying Lemma [61(b) yields 

conv (Pi U P_i) P P = {^Ri + ^ {*^} - -^0, 


where the inclusion follows from Pi -|- P_i C 2Pq 
one has relint(P P P) P Z'^ = relint(Po) P Z'^. By (a) 
have relint (Pq) P Z*^ = 


Hence, P P P = Pq. Therefore, 
Pq is lattice-free and hence, we 
X {—1,0,1}, where P P x {—1,1} 


Since P P Z'^ C Z-^-^ 

is contained in the boundary of P, this shows that P is a lattice-free integral polytope 
containing P. Since P is Z'^-maximal, we have P = P and in particular, Pi = Pi and 
P_i = P_i. 

Assertion (b3): Let v G R*^”^ be a non-integral vertex of Pq and let v := (u, 0) G R'^. 
Then v is not integral and since P is an integral polytope, v 0 vert(P). On the other 
hand, since u is a vertex of Pq, T cannot be written as convex combination of any points 
of PPP. Thus, we have v ^ A and consequently, in view of Q, we also have v G PPP. 
By m, there exist xi G Pi, x_i G P_i such that v = ^xi -|- ^X-i. Accordingly, 
2i; G Pi -|- P_i. This proves the assertion. □ 
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Analysis of special cross-sections 

From Figured on page [TOl one can make the following observation: three of the four 
maximal half-integral polytopes are simplices with the property that all of their integral 
vertices are unimodular. It turns out that if Pq in the setting of Theorem [7] is a simplex 
with this property, one can strengthen the assertions about the relation between Pi, P_i 
and Pq presented in Theorem [71 

Lemma 8 . Let Q G V{'L^) he a polytope such that Q is not a singleton and let v be a vertex 
of Q. Let g := gcd ( {(u), Cj) : w G veit{Q — v), i G {1,... ,d}} ), i.e., g is the greatest 
common divisor of all the components of the vertices ofQ—v. LetQ := ^{Q—v) G V{U^). 
Then every non-negative homothetic copy Q' of Q belonging to has the form 

Q' = aQ + b, where a G N U {0} and b G 

Proof. Let Q' G V{Z^) be a non-negative homothetic copy of Q. Since obviously Q' 
is then also a non-negative homothetic copy of Q, there exist a > 0 and 5 G such 
that Q' = aQ + b. By dehnition of Q, one of its vertices is o. Hence, b G yevt{Q') and 
therefore, in view of Q' G we have b G Z'^. Thus, Q' — b G P(Z'^) and hence from 

vert(Q' — b) = |( vert((5) — v)G'L'^ we get | {w' — v , ef) gZ for every w' G vert((5) and 
every i G {1, ... ,d}. Suppose now that a is not an integer, i.e. a has a representation 
as p/q with p G N and q G N \ {1} being relatively prime. Then q divides r,ei) 
for every w' G vert((5) and every i G {1,... ,d} and hence, qg divides {w , e*) for every 
w G vert((5 — v) and every i G {1, • • • ,d}, which contradicts the definition of g. □ 

Theorem 9. Let P, Pq, Pi, P_i be as in Theorem\^ If Po 'Is a (d—l)-dimensional simplex 
such that all integral vertices of Pq are unimodular, then the following statements hold: 

(a) Pi and P-i are non-negative homothetic copies of Pq, 

(b) Pi + P_i = 2 Po. 

Proof. Assertion (a). Let ui,... ,Ud G Z‘^“^\{o} be primitive vectors such that U{Pq) = 
{mi, ... ,Ud}. We circumscribe a non-negative homothetic copy of Pq around Pi. More 
precisely, we consider the polytope 

P{ := |x G : {x , Uj) < h{Pi,Uj) for every j G {1, ... , d}| , 

which contains Pi and is a non-negative homothetic copy of Pq since P[ is a simplex 
having the same set of outer normal facet vectors as Pq. We will now prove that P[ = Pi, 
which implies that Pi is a non-negative homothetic copy of Pq. 

Let V G vert(P{). Appropriately reindexing ui,... ,Ud, we assume 

{v , Ui) = h{P{,Ui) for every i G {1,... ,d — 1}. 

Denote by vq the vertex of Pq with the property 

(uo , Ui) = h{Po,Ui) for every i G {I,... ,d - 1]. 

Let us hrst consider the case vq 0 Z^~^. In this case, (|7|) yields 2vq G Pi -\- P_i. Thus, 
taking into account that 2vo is a vertex of 2Pq and Pi -|- P_i C 2Po, we get that 2vo 
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is a vertex of Pi + P_i. Note that Pi + P_i = conv(vert(Pi)) + conv(vert(P_i)) = 
conv(vert(Pi) + vert(P_i)) (see |Sch93t Theorem 1.1.2]) and hence, vert(Pi + P_i) C 
vert(Pi) + vert(P_i). Therefore, there exist vi G vert(Pi), V-i G vert(P_i) with 2vo = 
vi + V-i. In view of Pi + P_i C 2Po, we have cone(wi,..., ri^-i) = N{2Po,2vo) C 
N{Pi +P_i,2i;o). Furthermore, we have N{Pi + P_i,2i;o) = N{Pi,vi) n A^(P_i, u_i); 
see |Sch93[ Theorem 2.2.1]. Together, this yields cone(Mi,..., Ud-i) ^ N{Pi,vi). In 
other words, ui,..., Ud-i are outer normal vectors of Pi at vi. Hence, we have 


{vi , Ui) = h{Pi,Ui) for every i G {1,..., d - 1}, 

and thus, v and vi are defined by the same (uniquely solvable) system of linear equalities. 
Hence, v = vi and in particular, v is integral. 

We now treat the case that uq G i.e., vq is a vertex of Pq. In this case, since all 

integral vertices of Pq are unimodular, ui,..., u^-i form a basis of Z‘^“^, i.e., the matrix 
with rows Mi,...,Urf_i has determinant one. Since the vertex v of P{ is the unique 
solution of the system of linear inequalities given by 


{ui, v) = h{Pi,Ui) for every i G {1,... , d - 1}, 

and h{Pi,Ui) is integral for every i G { 1 , ..., d — 1 }, it thus follows that v is integral. 

We have now shown that all vertices of P{ are integral and hence, P{ is an integral 
polytope containing Pi. Observe now that P[ +P_i C 2Po, since for x G P( and y G P_i 
and every j G { 1 ,... , d}, one has 


{x + y, Uj) = {x, Uj) + {y, uj) < h{Pi,Uj) + h{P_i,Uj) = h{Pi + P-i,Uj) < h{ 2 Po,Uj). 


Here, the last inequality follows from Pi + P_i C 2Po. In view of Theorem [ ^b2) with 
Pi = P{ and P_i = P_i, we then immediately obtain Pi = P{- This proves that Pi is 
a non-negative homothetic copy of Pq. To show that P_i is a non-negative homothetic 
copy of Pq, one can argue in exactly the same way by interchanging the roles of Pi and 
^- 1 - _ 

Assertion (b). Let v be any vertex of 2Po. We define Pq := |(2Po — v) using g as 
in Lemma [H] with Q = 2Pq. By (a), Pi and P_i are homothetic copies of Pq. Thus, 
Lemma E] yields Pi = aPo -|- b and P_i = o'Pq + 6 ', where a, a' G NU {0} and b, b' G 
By Theorem l ^bl) we have Pi -|- P_i C 2Po. Summarizing these arguments, we get 


aPo -|- o!P q -|- (6 -|- h') C 2Po — gPo + v. 

Since g, a, a' > 0 and since Pq is (d—l)-dimensional, we have a-|-a' < g, i.e. 0 < o' < g—a. 
Then since 5, a G NU {0}, the polytope P_i := {g — a)Po -|- b' is in P(Z'^“^) and contains 
P_i. In view of Pi -|- P_i C 2Po and Theorem [ ^b2) we obtain P_i = P_i and by this, 
a' = g — a and a + a' = g. This yields gP^ -|- (6 -|- h') = gP^ + v and implies b + b' = v. □ 


Enumeration for lattice width two and dimension three 

In view of the results of this section, all possible choices for Pq are contained in the list 
of polygons given in Theorem [H Based on this, we can now construct all Z^-maximal 
lattice-free polytopes P in P(Z^) with lattice width two. We determine all such P and 
see that all of them are also R^-maximal. 
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Proposition 10. Let P G V{’l?) be a lattice-free polytope with lw(P) = 2. Then, P is 
-maximal if and only if P is unimodularly equivalent to one of the following polytopes 
(see Figure 

the simplex := conv(— 2 ei,4ei,ei + 862 , 263 ), 
the simplex := conv(o, 46i, 462 , 263 ), 

the simplex '■= conv (—61 + 62,61 + 862 , 263 , 26 i — 262 + 263 ), 

the simplex M 44 := conv(— 61 , 861,61 + 462 , 263 ), 

the pyramid M 5^4 := conv( 6 i - 62,-61 + 62 , 861 + 62,61 + 862 , 263 ), 

the prism '■= conv(- 6 i, 61 , 262 , 263 , 26 i + 263 , 61 + 262 + 263 ), 

the parallelepiped Mg ,2 := conv(o, —61 + 62 , 262,61 + 62 , 263,61 + 62 + 263,261 + 
263,61 - 62 + 263). 

In particular, P is -maximal if and only if P is -maximal. 

Proof. We remark that throughout the proof, whenever we claim that two polytopes 
which are explicitly given are unimodularly equivalent, we omit the verification of this 
claim. This verification can be carried out either by hand, using elementary linear 
algebra, or using specialized computer software. 

Let P be Z^-maximal. By applying a unimodular transformation, we can assume 
P C X [—1,1] and w(P, 63 ) = 2. For i G { — 1,0,1}, let Pi := {x G : (x, i) G P}. 
In view of Theorem 0] and Theorem [ ^a)[ we have that Pq is unimodularly equivalent to 
one of the polygons Q 2 ,Q 3 ,Qi,Q 5 as in Theorem 01 By applying another unimodular 
transformation to P, we assume that Pq = Qi for some i G {2,8,4, 5}. Furthermore, Pi 
and P_i satisfy the conditions of Theorem [ ^b)[ 

Case 1: Pq is a simplex. That is, P is Q 2 , Q 3 or Q 4 . These latter three triangles have 
the property that each of their integral vertices is unimodular. In view of Theorem [9l 
we obtain that if Pq G {Q 2 ,Q 3 ,Q 4 }, then Pi,P_i G P(Z^) are non-negative homothetic 
copies of Pq such that Pi + P_i = 2Po. In order to apply Lemma [H we fix a vertex 
V G P(Z^) of 2Po and consider Pq := ^(2Po — v), where g := gcd({t(;i, 1 x 2 , 163 , 1 x 4 }) 
with {wi,W 2 ) and ( 163 , 1 x 4 ) being the elements of vert(2Po — x) \ {o}. Then there exist 
a, o' G N U {0} and b, b' G Z^ such that Pi = oPq + 6 , P_i = o'Pq + b' and 2Po = 
gPQ + X = (o + a')Po + (6 + b'). As shown in the proof of Theorem [9](b), we have 
b + b' = V, i.e. b' = V — b, and a -\- a' = g. It is also obvious that reversing the roles of 
Pi and P_i corresponds to reflecting P with respect to x {0}, which is a unimodular 
transformation. Without loss of generality, we can therefore assume o < o'. Moreover, 
observe that we can also assume b = o (and hence, 6 ' = x). To see this, let b = {bi, 62 ) and 
let —>■ be the linear mapping given by (/?(xi, X 2 , X 3 ) = (xi — 61 X 3 , X 2 — 62 X 3 , X 3 ). 

Then ip is indeed a unimodular transformation mapping (aPo + 5) x {1} onto aPo x {1} 
and (a'Po + x — 6 ) x {—1} onto (o'Pq + x) x {—1} while leaving qPq x {0} unchanged. 

We now proceed as follows; for given Pq, we select a vertex x and determine g and 
Pq. Then for every pair of non-negative integers a, a' with a < a' and o + a' = g', we 
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check the polytope P = conv ((a'Po + f) x {—1} U oPq x {1}) for Z^-maximality. If P 
is Z^-maximal, we show that it is unimodularly equivalent to one of the seven polytopes 
given in the formulation of the theorem. 

For the three cases, we select the vertex v and obtain Pq as follows: 


Case 1.1: 

To = Q2, 

II 

0 

Po = conv ((0,0), (1,0), (0,1)) 

Case 1 . 2 : 

Pq = Qs, 

u = (- 2 , 0 ), 

= conv (( 0 , 0 ), ( 2 , 0 ), ( 1 , 1 )) 

Case 1.3: 

II 

u = (- 1 , 0 ), 

= conv (( 0 , 0 ), ( 2 , 0 ), ( 1 , 2 )) 





Case 1.1 Case 1.2 Case 1.3 

Figure 6 : Illustration of Pq and Pq (shaded) in the Cases 1.1-1.3 of the proof of Propo¬ 
sition [ini 


Case 1.1: Pq = Q 2 - Then 2Po = 4Po. In view of a < a', we have three cases to dis¬ 
tinguish: first, a = 0 and a' = 4, second, a = 1 and a' = 3, and third, a = a' = 2. In the 
first case, P is the simplex M 4^4 — 63 . In the second case, P is conv ((3Po —e 3 )U(Po-|-e 3 )), 
which is properly contained in the integral lattice-free polytope conv(o, 3ei, 3e2, 363 ) — 63 
and is thus not Z^-maximal. If a = 2, then P = conv(o, 2ei, 2 e 2 ) x [—1,1], which is 
properly contained in the integral lattice-free polyhedron conv(o, 2 ei, 2 e 2 ) x M and thus 
not Z^-maximal. 

Case 1.2: Pq = Qq. Then 2Po = 3Po — (2,0). There are two cases to distinguish: 
a = 0, a' = 3 and a = 1, a' = 2. In the former case, P is the simplex ~ ^ 3 . In the 
latter case, 

P = conv ((-2, 0, -1), (0, 0, -1), (-1,1, -1), (0,0,1), (2,0,1), (1,1,1)), 

which is properly contained in conv ((—2,0, —1), (0,0, —1), (—1,1, —1), (0, 0, 2)). This 
implies that P is not Z^-maximal. 

Case 1.3: Pq = Q^. Then 2Po = 2Po — (1,0). There are two cases to distinguish: 
o = 0, a' = 2 and a = a' = 1. In the former case, P is the simplex M 44 — 63 . If 
o = a' = 1, then P = conv ((Pq — (1, 0)) x {—1} U Pq x {1}) = — 63 . 

Case 2: Pq = Q 5 . Then we have 2Po = conv(±(2,0), ±(0, 2)) -|- (1,1). Since Pq 
does not have integral vertices, combining ([6]) and © we obtain Pi -|- P_i = 2Po. 
Consequently, for every u € Z^ \ {o}, one has F{2Pq,u) = P(Pi -|- P_i,m) = P(Pi,u) -|- 
P(P_i,tt); see |Sch93) Theorem 1.7.5(c)]. The latter equality shows that every one¬ 
dimensional face of Pi and P_i is parallel to an edge of Pq. In other words, if Pi or P_i is 
not a singleton, each of its facets is parallel to one of the line segments 5i := [( 0 , 0 ), ( 1 , 1 )] 
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and S 2 := [( 0 , 0 ), ( 1 , — 1 )] and thus, there exist non-negative integers ki,k 2 ,li,l 2 and 
vectors v^w ^ 1 ? such that 

Pi = kiSi + k 2 S 2 + V and P-i = /i5i -|- I 2 S 2 + w. 

Observe also that one has 2 Pq = 25i -|- 2 S 2 + (—1, 1), which in view of Pi -|- P_i = 2Po 
implies v + w = (—1, 1) and ki + li = k 2 + I 2 = 2. As argued in Case 1, reversing the 
roles of Pi and P_i results in reflecting P with respect to x {0} and hence, we can 
assume ki < li, which in view of fci -|- ^i = 2 implies ki < 1 . Likewise, in Case 1 we have 
fixed b = (0,0) and here, we can do the same for v, which implies w = (—1,1). Since li 
and I 2 are determined by the choice of ki and k 2 , respectively, we have the following six 
cases depending on the choice of (fci, /C 2 ) G {0,1} x {0,1,2}. 

• If (fci, /C 2 ) = (0, 0), then P = conv(2Po x {—1} U (esj) = — 63 . 

• If (fei, /C 2 ) = (0,1), then P = conv ((2Si -|- ^2 -|- (—1,1)) x {—1} U ^2 x {1}). Here, 
P is a prism unimodularly equivalent to M^^ 2 - 

• If (fci, /C 2 ) = (0, 2), then P = conv ((25i-|-(—1,1)) x {- 1 }U( 252 ) x (Ij) = ^ 4 ^ 2 - 63 - 

• If {ki, ^ 2 ) = (1,0), then P is a prism unimodularly equivalent to the case (fci, /C 2 ) = 
( 0 ,1) and hence to M^^ 2 - 

• If {ki,k 2 ) = (1,1), then P = conv ((Si -I-S 2 -|- (-1,1)) x {-1} U (Si -|- S 2 ) x {1}) = 
Ale ,2 — 63 . 

• If {ki, ^ 2 ) = (1, 2), then P is a prism unimodularly equivalent to the case (fci, /C 2 ) = 
( 0 ,1) and hence to M^^ 2 - 

This completes the proof that if P is Z^-maximal, it is unimodularly equivalent to 
one of the seven polytopes specified in the proposition. Conversely, the seven polytopes 
listed are all M^-maximal (since each of their facets is blocked) and, in particular, also 
Z^-maximal. □ 

4 The case of lattice width at least three 

In view of Proposition fTOl to complete the proof of Theorem [T] it remains to show that 
every Z^-maximal lattice-free polytope P G P(Z^) with lattice width at least three is also 
M^-maximal. Our goal is to derive an algorithm for a computer search which finds all 
such polytopes and verifies the latter claim. To be able to do so, we first need reasonably 
good bounds on the (suitably defined) size of P. 

Bounds on the size of lattice-free polytopes in P(Z^) with lattice width 
at least three 

The following parameters will be used to quantify the size of a lattice-free polytope 
P G P(Z^). We will use the volume of both P and its difference body P — P and also 
the lattice diameter of P. Finally, we will also use the so-called successive minima: for a 
convex body C C centrally symmetric with respect to the origin, i G {1,... , dj, the 
minimal value A > 0 such that XC contains i linearly independent points of Z'^ \ { 0 } is 
denoted by Xi{C) and called the i-th successive minimum of C. We remark that clearly, 
0 < Ai(C') < ... < Xd{C) and also Ai(tC') = lAi(C') for every positive number t. We 
will quantify the size of P in terms of 
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In order to bound the volume from above, we make use of the first and second 
theorem of Minkowski; see [Gru931 Theorems 4 and 12], which link the volume of an 
o-symmetric body to its successive minima. For a convex body C C symmetric with 
respect to the origin, one has 


vol(C') < 


< 


ntiAi(c’) “ ^i(c) 


d' 


( 11 ) 


Consequently, lower bounds on the first successive minimum imply upper bounds on the 
volume. Let iL be a convex body in not necessarily a centrally symmetric one. The 
volume of K and the volume of its difference body K — K are related by the well-known 
inequality 2“^ vol(iL) < \o\{K — K)] see for example |Sch93[ Theorem 7.3.1]. Since K — K 
is o-symmetric, one can apply m to obtain 


vol(iL) < 


1 


1 

- Xi{K -Ky 


( 12 ) 


In the following, we will derive a lower bound on Ai(P — P) for a lattice-free polytope 
P with lw(P) > 3. 

To establish a connection between the lattice width of a convex body and the suc¬ 
cessive minima of its difference body, we use another sequence of parameters connected 
to a convex body K. For i E {1,... , d}, we denote by the i-th covering minimum 

of K, defined to be the infimum over all p > 0 for which + Z'^ meets every {d — i)- 

dimensional affine subspace of (see [KL88j ). Clearly, 0 < ldi{K) < ■ ■ ■ < iJ,d{K). 
Kannan and Lovasz |KL88l Lemmas 2.3 and 2.5] provided a link between covering min¬ 
ima and successive minima by showing that 


+ Xd-i{K - K) (13) 

holds for every i E {l,...,d — 1}. Furthermore, they observed Iw(iL) = 

Note that the value fJ-d{P) is the standard covering radius (also called inhomogeneous 
minimum) of K and that UdiK) > 1 if and only if some translation of K is lattice-free. 

Hurkens [HurhOj proved the following relation between 1^2 and for iL C in the 
case d = 2: 

/i2(P)< (l + ^)w(i^)- ( 14 ) 

Kannan and Lovasz [KL88j observed that this also implies that (11411 holds for every 
d >2. Let us now return to the case d = 3. For a lattice-free convex body K C we 
can therefore combine (m with (|13p for i = 2 and the fact that fJ- 3 {K) > 1 to obtain 
the following chain of inequalities: 

^ll{K) + Xl{K -K). 

Consequently, one has 

Ai(P-K)>l-(l + ^)/ri(iL). (15) 


1 < hsiK) < /i2(P) + Xi{K -K)<{1 + 


V3 


21 


















Recall that is the reciprocal of Iw(Rr). Hence, the right hand side of (1151) is 

positive whenever Iw(it') > 1 + Applying these observations to the case of a lattice- 
free polytope P G P(Z^) with lattice width at least three we obtain a positive lower 
bound on Ai(P — P) from (fT^ : 


We are now ready to list the inequalities needed for the proof of Theorem [TJ 

Proposition 11 (Size bounds for lattice-free polytopes in P(Z^) with lattice width 
at least three). Let P G P(Z^) he a lattice-free polytope with lw(P) > 3. Then, the 
following statements hold: 

(a) vol(P) < 27, 

(h) vol(P-P) <8-27, 

(c) Ai(P-P)>i, 

(d) ld(P) < 3. 

Proof of Proposition 1771 Let Aj := Ai(P —P) and := fii{P) for i G {1,2,3}. Assertion 

(c) follows immediately from (|16p . 

We now prove (a) and (b). By (flTT) and (fT^ . it suffices to prove that the product 
A 1 A 2 A 3 is bounded from below by 1/27. In the case Ai > 1/3, this is clear, since in this 
case A 1 A 2 A 3 > Af > 1/27. We switch to the case Ai < 1/3. Applying (fT^ for i = 1, 
we get X 2 > IJ -2 — hi applying (fT3l) for i = 2 yields h2 P h3 ~ ^i- Combining these 
two inequalities, we get A 2 > ps — Xi — pi. Since P is lattice-free, we have p^ > 1. 
We further have pi < 1/3 because l/pi = lw(P) > 3. Hence, A 2 > 2/3 — Ai and 
thus, A 1 A 2 A 3 > Ai(2/3 — Ai)^. Taking into account (c), we have 1/4 < Ai < 1/3. 
For a moment, we view Ai as a variable and determine a lower bound of the function 
/(Ai) := Ai(2/3 — Ai)^ on the interval (|,|]- The latter function is decreasing on 
this interval and hence attains its minimum on this interval for Ai = 1/3. This yields 
A 1 A 2 A 3 > /(Ai) > /(1/3) = 1/27 and completes the proof of (a) and (b). 

This leaves us to show (d). We observe the following connection between lattice 
diameter and first successive minimum. By the definition of ld(P), there exist distinct 
points 2 , 2 ;' G P n with u := (-z — z') G Z^. Thus, u is a non-zero integral vector 

in id(p) {P ~ P)- Taking into account the definition of the first successive minimum, we 
get 

Hence, we have ld(P) < 1/Ai < 4 in view of (c) and since ld(P) G N, this shows 

(d) . □ 
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Completing the case of lattice width at least three by a computer search 

Based on the bounds on the size of lattice-free polytopes in P(Z^) with lattice width at 
least three, we can now derive a search algorithm that finds, up to affine unimodular 
transformations, all polytopes in the set 

{P G V{'I?) : P is Z^-maximal lattice-free and lw(P) > 3} . (17) 

We present the algorithm along with mathematical arguments that verify its correctness. 
We give a rather high-level description omitting algorithmic details that can be resolved 
in a straightforward manner. We implemented our search algorithm in Python 2.7 - 
a high-level language with a self-explanatory syntax. The code of this implementation 
can be found in the appendix. Our implementation contains a straightforward test that 
checks for each polytope P in (jl7l) whether it is also M^-maximal; we check whether the 
relative interior of each facet of P contains an integral point. As a result, an execution 
of our code confirms (within less than half an hour) that every polytope in (I17p is indeed 
M^-maximal. 

Since all M^-maximal lattice-free polytopes in P(Z^) were enumerated in |AWWli] . 
it was not our intention to output an irredundant list of all polytopes in (I17p . In fact, 
our algorithm enumerates many polytopes that are pairwise unimodularly equivalent. 
We decided to refrain from excluding such redundancy to keep the code as simple as 
possible and thus more easily verifiable. 

Step 1: Fixing a boundary fragment 

Let P € P(Z^) be Z^-maximal lattice-free and let i := ld(P). In view of Proposition [Til 
the possible values of I are 1, 2, and 3. Because P is lattice-free, there exists a facet F of P 
with ld(F) = ld(P) = I and by applying an appropriate unimodular transformation, 
we can assume F P M? x {0}. Let us first treat the case that I = \. Howe’s Theorem 
(see, e.g., |Sca85[ Theorem 1.5]) asserts that if P G P(Z^) satisfies P r\T? = vert(P), 
then the lattice width of P is one. Since, up to unimodular equivalence, x [0,1] is 
the only Z^-maximal lattice-free polyhedron with lattice width one, such P is not Z^- 
maximal. Hence, P has to contain an integral point which is not a vertex of P, and 
since P is lattice-free, this point has to be in the boundary of P. Clearly, there cannot 
be an integral point in the relative interior of an edge of P, since otherwise this edge 
has at least three collinear integral points, contradicting I = \. This implies that we 
can choose P such that relint(P) n Z^ / 0. Suppose |P n Z^j > 4. Then two points 
in |P n Z^l, say X and y, are congruent modulo two and hence, the line segment [x, y] 
contains ^(x -|- y) G Z^, a contradiction to = 1. Thus, P is a triangle with precisely 
one integral point in its relative interior. It follows that P is unimodularly equivalent 
to conv((-l,-l,0), (1,0,0), (0,1,0)) ; see |R,ab89j . where all integral polygons with 
precisely one interior integral point were listed. 

It is easy to check that in the the case £ = 2, up to unimodular transformation, P 
contains the triangle conv ((0,0,0), (2,0,0), (0,1, 0)) and likewise, in the case f = 3, up 
to unimodular transformation, P contains the triangle conv ((0,0,0), (3,0,0), (0,1,0)). 
Thus, we may assume that P contains one of the following three triangles B in its 
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boundary: 


B := 


'conv((-l,-l,0), (1,0,0), (0,1,0)) 
< conv ((0,0,0), (2,0,0), (0,1,0)) 
,conv ((0,0,0), (3,0,0), (0,1,0)) 


if ^ = 1, 
if ^ = 2, 
if ^ = 3. 


(18) 


Step 2: Inscribing a pyramid 

Clearly, after fixing (. and the corresponding B, without loss of generality we can assume 
that we have X 3 > 0 for all points {xi,X 2 -,x^) E P. Let h := max{x 3 : (xi,X 2 ,a: 3 ) E P} 
and fix some vertex of P of the form a = {ai,a 2 ,h). By the definition of the lattice 
width, we have h > lw(P) and since we assume lw(P) > 3, we have h > 3. By 
applying an appropriate unimodular transformation that keeps the third component 
unchanged, we may assume that 0 < ai, 02 < /i — 1 holds. To see that, observe that for 
every pair of integers k,k' the linear mapping if given by ip{ei) = ei, ip{e 2 ) = 62 and 
= kei + k'e 2 + 63 is unimodular and maps (ai, 02 , h) to (ai + kh, 02 + k'h, h). In 
order to obtain an upper bound on h, we define the tetrahedron T := conv(B U {a}) and 
observe that we have 

(®) • vol(r) = vol(r - T) < vol(P -P)<8 - 27, 

where the first eqnality follows from |RS57l Theorem 2] and the last inequality fol¬ 
lows from Proposition [TT] Denoting by area(B) the area of B and using vol(T) = 
h- area(B), we thus obtain 

, 3-8-27 

h < - • 

( 3 )- aiea{B) 

In the cases i = 2 and £ = 3, this bound evaluates to /i < 32 and h < 21, respectively. 
However, in the case i = 1, the following Lemma provides a tighter upper bound on h. 

Lemma 12. Let P E R(Z^) he a lattice-free polytope such that P contains (oi, 02 , h) E 
and also the points (—1,—1,0), (1,0,0) and (0,1,0). Then, h< 12. 

Proof. We argue by contradiction by assuming L > 12 and showing that this contradicts 
the lattice-freeness of P. Let D denote the triangle conv ((—1, —1, 0), (1, 0,0), (0,1, 0)) 
and consider D' := D Ci {—D). Observe that o & D and hence, o E D'. Furthermore, D' 
is symmetric with respect to o by construction. We also introduce the double pyramid 
R := conv(Zl' U {±(ai, 02 ,/i)}). One can verify directly that the area of D' is one and 
thus, the volume of R is |/i > 8. Since R is an o-symmetric convex body, we can apply 
IfTTI) to obtain > vol(i?) > 2^ or, equivalently, Xi{R) < 1. By definition of the 

snccessive minima, this implies that int(i2) contains a point z of \ {o}. Because of 
the symmetry of R, this also implies —z E int(i7). As D' n Z^ \ { 0 } = 0, either z or —z 
has to be contained in int(i?) O (M^ x (0, 00 )). The latter set, however, is contained in 
int(P), a contradiction to the lattice-freeness of P. □ 

We now summarize the previous observations as follows: we may assume that P is 
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contained in x [0,/i], where 


h : = 


'12 
< 21 


if £ = 1, 
if £ = 2, 


32 if ^ = 3, 

\ ■ 


and that P contains the tetrahedron T = conv(il U {a}), where B is defined as in (|18l) 
and 


a ^ := {(xi,X 2 ,X 3 ) E : 0 < xi,X 2 < X 3 — 1, 0 < X 3 < h} . (19) 

We now enumerate all elements of the set in (I19p . Without further qualifications, this 
enumeration yields several ten thousand points. However, recall that P is lattice-free and 
satisfies Xi{P — P) > j by Proposition [TTl Since TCP, these two properties also have 
to be satisfied by T. Furthermore, because P has lattice diameter I and T C P, we have 
ld(T) < i, while on the other hand B CT yields ld(T) > Id(il) = i. Hence, we obtain 
that T is lattice-free and satisfies ld(T) = I as well Ai(T —T) > The computer search 
checks for each of the possible choices for a = ( 01 , 02 , h) whether these three properties 
are fulfilled for the tetrahedron conv(H U {(oi, 02 ,/i)}). As a consequence, we end up 
with only 69 valid choices for a in total. 

Step 3: Collecting all candidate vertices 

In the next step, for a given combination of fixed lattice diameter i and pyramid T 
(given by B and o = ( 01 , 02 ,/i) as introduced above), we enumerate all lattice-free 
polytopes P in P(Z^) which have lattice diameter £ and contain T. In order to do so, we 
enumerate a finite subset of Z^, depending only on T, which we denote by cand(T) (we 
call the elements of cand(T) candidate vertices) such that for every P as above one has 
vert(P) C cand(T). This is done as follows. First, recall that every vertex v = {vi,V 2 , 03 ) 
of P satisfies 0 < U 3 < by our choice of h. By Proposition 1111 we further have 

vol (conv(r U {u})) < vol(P) < 27. 

The following lemma shows that all points x E satisfying vol (conv(T U {x})) < 27 
lie in an appropriate homothetic copy of T. 

Lemma 13. Let T C be a tetrahedron, c the barycenter of its vertices and let t > 
vol(r). Then 

|x E : vol ( conv(T U {x})) < t| C A- T -|- (1 — A)- c, 
where X := 4- +1. 

Proof. Consider a point x E not belonging to XT -|- (1 — A)c. We show that then one 
has vol (conv(T U {x})) > t. Let /3i,... ,/34 be the barycentric coordinates of x with 
respect to the vertices of XT -|-(1 —A)c. It is straightforward to verify that the barycentric 
coordinates of x with respect to the vertices of T are /?i -|- ...,[34 + Since x 0 
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AT + (1 — A)c, at least one of /3i,..., (34, say /3i, is negative and hence, the corresponding 
barycentric coordinate of x with respect to T is 71 := /3i + + 1 < 0. 

Let vi be the vertex of T associated with this barycentric coordinate and let T ;= 
conv(vert(T) \ {fi}). Let denote the distance of vi from H := aff(vert(T) \ {wi}) 
and let dx denote the distance of x from H. Then from the geometric interpretation of 
the barycentric coordinates, one has dx = | 7 i|di;i = —'Jidvi- Likewise, it follows from 
the geometric interpretation and the fact that 71 < 0 that T and conv(T U {x}) do not 
have common interior points. In view of —71 < 1 — , this yields the following chain 

of inequalities: 

vol (conv(T U {x})) > vol(T) + vol ( conv(T U {x})) 

= vol(T) — 7 i vol (conv(T U {xi})) 

= (1 - 71 ) vol(T) > vol(T) (1 + - 1 ) =t, 

which proves the statement. □ 

Thus, if c denotes the barycenter of the vertices of T, each vertex of P is contained 
in 

s(T) := {(xi,X 2 ,X 3 ) G (A-T + (1 - A)-c) n : 0 < X 3 < h} , 

where A := 4- ^ + 1- We call s(T) the search region for T. Similarly to the 

previous step, for each vertex x of P we have that conv(T U {x}) is lattice-free and its 
lattice diameter is bounded by £. Hence, a valid choice for the set cand(T) is 

cand(T) := {x G s(T) : conv(T U {x}) is lattice free, 

ld(conv(T U {x})) = f\. 

For the enumeration of cand(T), we enumerate s(T) and check the defining conditions 
for cand(T). Since for most choices of T, the set s(T) is huge, this is the most time- 
consuming step. Therefore, in order to reduce the running time, we define cand(T) 
without use of the successive minima. At the end of Step 3, cand(T) is stored as a list 
of points and available for further computations. 

Step 4: Combining the candidate vertices 

In the final step, for each fixed £ and T it remains to enumerate all lattice-free polytopes 
with lattice diameter I and lattice width at least three which can be obtained by taking 
the convex hull of T and a set of elements of cand(T), i.e., all polytopes in the set 

S:={P £ V{Z^) ; vert(P) C cand(T), T C P, P lattice-free, ld(P) = £, lw(P) > 3} . 

This is done in the following way. Let xi,..., x^ be the points of cand(T). We construct 
S iteratively by listing the elements of sets Sq, ..., Sk, where Sq = {T} and Sk = S. The 
set Si is constructed from 5j_i in a way that for every i G {0,..., A:}, 5* is the set of all 
polytopes in P(Z^) satisfying the following conditions; 

vert(P) C vert(T) U {xi,... ,Xi}, 

T C P, P is lattice-free, ld(P) = £, 
lw(PU {xi+i,... ,Xfc}) > 3. 
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( 20 ) 

( 21 ) 

( 22 ) 



Obviously one has Sq := {T} and can compute Si as the union of the set of all polytopes 
of Si-i satisfying (| 2 ^ and the set of all polytopes P = conv(P'u{ui}), where P' G Si-i, 
satisfying conditions (ini) and (|22p . Only the third condition needs explaining: this 
is to ensure that for our preliminary polytope ‘built’ from T and the first i elements 
of cand(r), which might not have lattice width three, it is still possible to obtain a 
polytope of lattice width at least three by ‘adding’ points of cand(T) which have not yet 
been considered. If this is not possible, there is no need to consider this polytope any 
further. Observe that we have S = Sk and that each Si can be recursively computed for 
z G {1,... ,/c}. 

With an arbitrary choice of the sequence the cardinalities of some of 

the sets Si,...,Sk can be so large that a computation becomes practically impossible. 
However, this can be avoided by sorting the Vj's decreasingly according to the absolute 
value of their second coordinates. The idea behind this order is that, for most T, the 
majority of points in cand(T) (including the vertices of T) have a second coordinate with 
small absolute value. Thus, in order to arrive at a polytope with lattice width at least 
three, at least one of the Uj’s appearing early in the sorted list has to belong to P. Due 
to condition (j22p , any polytope not containing one of the first Vj ’s is then discarded from 
the sets Si for already small i. Furthermore, for any computed polytope containing one 
of the Uj’s from the beginning of the sorted list, the number of further candidates that 
can be added to arrive at a lattice-free polytope turns out to be very small. 

In the described procedure, for all constructed polytopes we need to decide whether 
their lattice width is at least three. This is done by simply computing the width with 
respect to a few explicit directions; more specifically, the directions in the set V defined 
in the following lemma. 

Lemma 14 . Let P be a polytope containing o, ei, 62 and some point a = (oi, 02, h) G 
with h > 1 . Then lw(P) > 3 holds if and only if 

w(P,v) >3 (23) 


holds for all v in the set 

V := {(ui,U2,U3) G Z^ \ { 0 } : |ui|, |t> 2 |, lui - U 2 I < 2, 

^(-2 - viai - 0202) <V 3 < i (2 - viai - ^202)}. 

Proof. If lw(P) > 3, then by definition of the lattice width, (I23p holds for every v G 
Z^\{o} and hence in particular for all u G H. To show the reverse implication, it suffices 
to show that if lw(P) < 2 holds, then there exists some u G H violating (1231) . Suppose 
we have lw(P) < 2. Then there exists a (primitive) vector v* = {vi,V2,v'^) G Z^ \ { 0 } 
with 


2> w{P,v*) > w({o, 61 , 62 }, u*) =max{|z;}|, bi-^^ 2 l}) 

where the second inequality follows from the fact that we have 0,61,62 G P. Finally, 
since we have o, a G P, we must have 

|(P,a) - (P,o)| < 2, 
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which is equivalent to 


— 2 < u^ai + ^202 + < 2, 

as claimed. □ 

Finally, having computed the set S^, for each polytope P G Sk we heuristically 
search for a certificate showing that P is not Z^-maximal. This is done by testing for a 
few points p G \ P in a certain neighborhood of P if conv(PU {p}) is still lattice-free. 
If this check fails for all considered choices of p, the polytope P is regarded as being 
potentially Z^-maximal and we then test whether P is even M^-maximal. It turns out 
that P is indeed M^-maximal in all these cases and hence, we obtain the following result. 

Proposition 15. Let P G P(Z^) be a lattice-free polytope with lw(P) > 3. Then, P is 
-maximal if and only if P is -maximal. 

Combining this with the results of the previous section, we obtain the proof of 
Theorem [H 

Proof of TheoremUl Let P G P(Z^) be lattice-free. If P is M^-maximal, then it is also 
Z^-maximal by definition of Z^-maximality. Conversely, if P is Z^-maximal and bounded, 
then in view of Propositions [101 and [15] we have that P is also M^-maximal. If P is Z^- 
maximal and unbounded, then it is unimodularly equivalent to conv(o, 2ei, 2 e 2 ) x M or 
to [0,1] X see [AWWlll Proposition 3.1]. In both cases, we have that P is also 
M^-maximal. □ 

Remark 16. Consider again the question asked in Remark [5l for which pairs s,fi G N 
are Z'^-maximality and M'^-maximality equivalent for lattice-free polyhedra in P(iZ'^)? 
Theorem [H shows that this equivalence does hold for d = 3 and s = 1. In Remark [5l 
we stated that equivalence also holds for d = 2 and s G {1,2}, while it trivially also 
holds for d = 1. For all other pairs of d, s G N, this equivalence does not hold: polytopes 
in P(iZ'^) which are Z'^-maximal but not R'^-maximal can be constructed following the 
ideas of [NZlll Proof of Theorem 3.2]. 
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Appendix 

This appendix contains an implementation of the algorithm presented in Section 0] in 
Python 2.7. Executing this code yields Proposition [T^ It is divided into three files, 
the first of which contains the main part of the algorithm, while the other two provide 
basic functions for computations on three-dimensional polytopes and basic functions and 
objects needed for geometric computations. 

search.py 


This is the main part of the code. The procedure is executed by calling the 
function "main". 


from basics import VectorSd 
from fractions import Fraction 
from polytopes import PolytopeSd 

import basics 
import math 

UPPER.BOUND.VOLUME = 27 

STRICT.LOWER.BOUND.SUCC.MIN = Fractiond, 4) 
def main O : 

"""Executes the search algorithm.""" 

# Step 1: enumerate boundary fragments (and run search) 
print "case: lattice diameter = 1" 

run.search(base = [Vector3d(x, y, 0) for (x, y) in [(1, 1), (—1, 0), (0, 
- 1 )]] , 

upper.bound.height = 12, 
lattice.diameter = 1) 

print "case: lattice diameter = 2" 

run.search(base = [Vector3d(x, y, 0) for (x, y) in [(0, 0), (2, 0), (0, 
1 )]] , 

upper.bound.height = 32, 
lattice.diameter = 2) 

print "case: lattice diameter = 3" 

run.search(base = [Vector3d(x, y, 0) for (x, y) in [(0, 0), (3, 0), (0, 
1 )]] , 

upper.bound.height = 21, 
lattice.diameter = 3) 
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33 # 


34 

35 def run_search(base, upper_bound_height, lattice.diameter): 

36 # Step 2: enumerate apices 

37 print " — computing apices . . . 

38 apices = list(enumerate_apices(base, upper_bound_beigbt, 

lattice_diameter)) 

39 print " — apices:", len(apices) 

40 

41 # Step 3: collect further candidate vertices 

42 for apex in apices: 

43 print " investigating pyramid with base", base, "and apex", apex 

44 print " — computing candidates ..." 

45 

46 pyramid = Polytopeld(base + [apex]) 

47 candidates = get_sorted_candidates(pyramid, lattice_diameter) 

48 print " — candidates:", len(candidates) 

49 

50 # Step 4: enumerate combinations 

51 #— searches for a ’critical’ polytope (see docstring of 

’search.critical.polytope ’ for a definition) with base 

52 # ’base’, apex ’apex’, and lattice diameter ’ lattice.diameter ’ 

53 # — as a result, no such critical polytope exists, which verifies 

the claim 

54 Q = searcb_critical_polytope(pyramid, apex, lattice_diameter, 

candidates) 

55 if Q: 

56 print "! error, found critical polytope", Q.vertices!) 

57 exit(l) 

58 

59 def enumerate_apices(base, upper_bound_beigbt, lattice_diameter): 

60 

61 Returns a list of all possible apices for the given boundary fragment 

’base’, see explanation of Step 2. 

62 """ 

63 for z in range(3, upper_bound_beigbt + 1): 

64 for X in range(z): 

65 for y in range(z): 

66 apex = Vector3d(x, y, z) 

67 pyramid = PolytopeSd(base + [apex]) 

68 

69 if is_valid_pyramid(pyramid, lattice_diameter): 

70 yield apex 

71 

72 def is_valid_pyramid(pyramid, lattice_diameter): 
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73 


74 Returns true iff ’pyramid’ is lattice—free, has lattice diameter of at 

most ’ lattice.diameter’ and has first 

75 successive minimum less than 1/4. 

76 

77 if not pyramid.is_lattice_free() : 

78 return False 

79 

80 if pyramid.has_ld_greater_than(lattice_diameter): 

81 return False 

82 

83 if pyramid . success!ve_minimum_at_most CSTRICT_LOWER_BOUND_SUCC_MIN): 

84 return False 

85 

86 return True 

87 

88 def get _sorted_candidates(pyramid, lattice_diameter): 

89 

90 Computes a sorted list of candidates for given pyramid ’pyramid’ and 

lattice diameter ’lattice.diameter ’ (see the 

91 explanations of Steps 3 and 4 as well as the remarks on the order of the 

list). 

92 """ 

93 candidates = [] 

94 

95 for c in integer_points_in_search_region(pyramid): 

96 if not pyramid.contains (c) : 

97 Q = PolytopeSd(pyramid.verticesO + [c]) 

98 if Q.is_lattice_freeO and not 

Q.has_ld_greater_than(lattice_diameter): 

99 candidates.append(c) 

100 

101 is_in_strip = lambda c: pyramid.ymin() <= c.y and c.y <= pyramid.ymax() 

102 

103 return [c for c in candidates if is_in_strip(c)] + [c for c in 

candidates if not is_in_strip(c)] 

104 

105 def integer_points_in_search_region(pyramid) : 

106 """ 

107 Computes a search region for integer points that can be added to 

’pyramid’ without exceeding the volume bound 

108 ’volume.bound ’ (see explanation of Step 3). 

109 """ 

110 vol = basics.volume_of_simplex(pyramid.vertices()) 

111 assert vol <= UPPER_BOUND_VOLUME 
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112 


113 

114 

115 

116 

117 

118 

119 

120 def 

121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 


vertices = pyramid.vertices() 

center = (vertices [0] + vertices [1] + vertices [2] + vertices [3]) * 
Fraction(1,4) 

scaling = 1 + 4 * (Fraction(UPPER_BOUND_VOLUME) / Fraction(vol) — 1) 
search_region = pyramid.homo.copy(scaling, center) 

return search_region.integer.points(zmin = 0, zmax = pyramid.zmax()) 

search_critical_polytope(pyramid, apex, lattice.diameter, candidates): 

This function computes all lattice—free polytopes with lattice diameter 
<= ’lattice^diameter ’ and lattice width >= 3 
that can be constructed from 'pyramid' by adding points of 'candidates'. 

For each such polytope, we check whether it 
is 'critical', i.e., potentially Z''3—maximal lattice—free but not 
R''3—maximal lattice—free. The method of generating 
such a list is described in the explanation of Step 4. 

if is.critical(pyramid, apex): 
return pyramid 

polytopes.built = set([pyramid]) 
num.candidates = len(candidates) 

while candidates: 

c = candidates.pop() 

polytopes.new = set(polytopes.built) 

# For each polytope P, keep P and add P' := P \cup \{ c \} if P' is 

lattice free and Id(P') <= lattice.diameter 
for P in polytopes_built: 
if not P.contains(c): 

Q = Polytope3d(P.vertices() + [c]) 
if Q.is_lattice_free() and not 

Q.has_ld_greater.than(lattice.diameter): 
if is.critical(Q, apex): 

return Q 

polytopes.new.add(Q) 

# Only keep polytopes that, by adding the remaining candidates, can 

be extended to ones with lattice width >= 3 
polytopes.built = set([P for P in polytopes.new if 


34 



149 


not 


has_lattice_width_at_most_two(P.verticesC) 
+ candidates, apex)]) 

150 print " — completed %d / %d (polytopes stored; %d)" % \ 

151 (num.candidates — len(candidates) , num.candidates, 

len(polytopes.built)) 


152 

153 

154 

155 

156 def 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 


if not polytopes.built: 

return 

is.critical(P, apex): 

Returns true iff P 

— has lattice width at least three, 

— is not R'^3—maximal lattice—free, but 

— seems to be Z'^S—maximal lattice—free. 

if has_lattice_width_at_most_twoCP■vertices(), apex): 
return False 

if P.all_facets_blocked(): 
return False 

# Heuristically try to add some integer points outside of P 
(itiin.x , min_y , itiin.z , max_x , max_y , max_z) = P. bounding_box () 
return not P.is_lattice_free_extendable([VectorSd(x, y, z) for x in 
range(min_x — 1, max_x + 2) \ 

for y in 

range (itiin_y 
- 1 , 
max_y + 

2 ) \ 
for z in 

range (itiin_z 
- 1 , 
max_z + 

2)]) 


173 

174 def has_l attice_width_at_most_two(points, apex): 

175 """ 

176 Returns true iff ’points’ has lattice width of at most two, see the 

lemma in explanation of Step 4. 

177 """ 

178 for vx in range(—2, 3): 

179 for vy in range(—2, 3) : 

180 if absCvx — vy) <= 2: 
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181 


min_z = basics.custom.ceil(Fraction(— 2 
apex.y , apex.z)) 

182 max_z = basics.custom.floor(Fraction (2 

apex.y , apex.z)) 

183 

184 for vz in range(min_z, max_z + 1): 

185 if vx ! = 0 or vy ! = 0 or vz ! = 0: 

186 V = Vector3d(vx, vy, vz) 

187 

188 

189 

190 return False 

191 

192 # - 

193 

194 mainO 

polytopes.py 

1 """Straightforward implementation of basic functions on three-dimensional 

(rational) polytopes.""" 

2 

3 import basics 

4 import itertools 

5 

6 from basics import VectorSd 
from fractions import Fraction 

8 

9 class Inequality3d(ob j ect) : 


10 

"""Represents a linear inequality of the form a * 

X <= b. """ 

11 




12 

def 

_init_(self, Ihs, rhs): 


13 


"""Creates a linear inequality of the form a * 

X <= b, where 



'Ihs' and b = 'rhs'.""" 


14 


self.Ihs = Ihs 


15 


self.rhs = rhs 


16 


if self.Ihs.is_zero(): 


17 


raise Exception("trivial inequality: %s" % 

str(self)) 

18 




19 

def 

__str__(self): 


20 


return str(self.Ihs) + " <= " + str(self.rhs) 


21 




22 

def 

_repr_(self): 


23 


return str(self) 


24 





if max([p * V for p in points]) — min([p * v for p 
in points]) <= 2: 
return True 


— vx * apex.x — vy * 

— vx * apex.X — vy * 
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25 


def __eq__(self, other): 

26 """Returns true iff all coefficients including the right-hand side 

of ’self’ and ’other’ coincide.""" 

27 return self.lhs == other.Ihs and self.rhs == other.rhs 

28 

29 def __ne__(self, other): 

30 return self.lhs != other.Ihs or self.rhs != other.rhs 

31 

32 def _hash_(self): 

33 return hash(str(self)) 

34 

35 def is_satisfied(self, p): 

36 """Returns true iff ’p’ satisfies the inequality.""" 

37 return self.lhs * p <= self.rhs 

38 

39 def is_ strictly.satisfied (self, p): 

40 """Returns true iff the inequality is strict for ’p’.""" 

41 return self.lhs * p < self.rhs 

42 

43 def holds_with_equality(self, p) : 

44 """Returns true iff ’p’ satisfies the inequality with equality.""" 

45 return self.lhs * p == self.rhs 

46 

47 # - 

48 

49 class Polytopeld(ohject): 

50 """Represents a three-dimensional (rational) polytope.""" 

51 

52 def _init_(self, points): 

53 """Constructs the convex hull of ’points’.""" 

54 assert(len(points) >= 3) 

55 

56 # Store bounding box 

57 self._xmin = min([p.x for p in points]) 

58 self._ymin = min([p.y for p in points]) 

59 self._zmin = min([p.z for p in points]) 

60 self._xmax = max([p.x for p in points]) 

61 self._ymax = max([p.y for p in points]) 

62 self._zmax = max([p.z for p in points]) 

63 

64 # Compute facet—defining inequalities 

65 #— this is done by iterating over all triples of points that span 

some plane H 

66 # — if all points in ’points ’ lie on one side of H (including H 

itself), this defines a facet 
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69 
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74 

75 

76 

77 

78 
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81 
82 

83 

84 

85 

86 
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88 

89 

90 

91 

92 
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97 

98 

99 

100 



self..facets = set() 

for triple in itertools.combinations(points, 3): 

normal = basics.norma1_vector_of .triangle(triple) 

if not normal.is.zero (): 

max.scal.prod = max([normal * p for p in points]) 
min.scal.prod = min([normal * p for p in points]) 

if max.scal.prod — min.scal.prod == Q: 

raise Exception("polytope low—dimensional") 

scal.prod = normal * triple [0] 

if max.scal.prod == scal.prod: 

self..facets.add(Inequality3d(normal, scal.prod)) 
elif min.scal.prod == scal.prod: 

self..facets.add(Inequality3d (—normal, —scal.prod)) 

# determine which points are vertices 

is.vertex = lambda p: len([f for f in self..facets if 
f.holds.with.equal!ty(p)]) >= 3 

self..vertices = frozenset([p for p in points if is.vertex(p)]) # 
frozenset important for hashing 

def bounding.box(self): 

"""Returns a tuple (xmin, ymin, zmin, xmax, ymax, zmax) containing 
the minimum and maximum coordinates of 'self 
with respect to all coordinate directions.""" 

return (self..xmin, self..ymin, self..zmin, self..xmax, self..ymax, 
self. .zmax) 

def integer.points.in.bounding.rect.of.slice(self, height): 

"""Yields a set S of integer points such that every integer point in 
'self with z—coordinate being equal to 
'height' is contained in S.""" 

# First, we compute the intersection of 'self with the set of 

points having z—coordinate being equal to 

# 'height'; note that every vertex of this set is 

# — a vertex of 'self with z—coordinate being equal to 'height', or 
#— the intersection of z = 'height' with a line segment between 

some vertex above and some vertex below z = 

# 'height ' 

above = [v for v in self..vertices if v.z > height] 
inside = [v for v in self..vertices if v.z == height] 
below = [v for v in self..vertices if v.z < height] 
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130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 


points = inside 

if len(above) > 0 and len(below) > 0: 
for a in above: 

for b in below: 

mu = FractionCa.z — height, a.z — b.z) 
points.append(a * (1 — mu) + b * mu) 

# By construction, ’self’ intersected with z = ’height’ is the 

convex hull of 'points'; now we compute a 

# bounding rectangle (xmin, ymin, xmax, ymax) of ’points’ 
if len(points) > 0: 

xmin = basics.custom.ceil(min([p.X for p in points])) 
ymin = basics.custom.cei1(min([p.y for p in points])) 
xmax = basics.custom.floor(max([p.X for p in points])) 
ymax = basics.custom.floor(max([p.y for p in points])) 

# Finally, we return all integer points inside the bounding 
rectangle 

for X in range(xmin, xmax + 1): 

for y in range(ymin, ymax + 1): 
yield Vector3d(x, y, height) 

def is_l attice.free(self): 

"""Returns true iff ’self’ does not contain any integer point in its 
interior .""" 

for height in range(basics.custom_ceil(self._zmin), 
basics . custoiti_floor (self . _zmax) + 1): 

for p in self.integer_points_in_bounding_rect_of_slice(height): 
if self.contains_in_interior(p): 
return False 

return True 

def contains_in_interior(self , p): 

"""Returns true iff ’p’ is contained in the interior of ’self’.’””’ 
for f in self..facets : 

if not f.is.strictly.satisfied(p): 
return False 
return True 

def contains(self , p) : 

"""Returns true iff ’p’ is contained in ’self’.""" 
for f in self..facets : 

if not f.is.satisfied(p) : 
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163 
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return False 
return True 

def ymin(self) : 

'""'Returns the minimum y—coordinate of points in the polytope 
’self ’. 

return sel£._ymin 
def ymax(self) : 

"""Returns the maximum y—coordinate of points in the polytope 
’self ’.""" 
return sel£._ymax 

def zmax(self) : 

"""Returns the maximum z—coordinate of points in the polytope 
'self ’.""" 
return self._zmax 


def integer.points(self, zmin = None, zmax = None): 

Returns all integer points of ’self’ with z—coordinate between 
’zmin’ (default: —inf) and ’zmax’ (default: inf). 

if zmin is None: 

zmin = basics.custom.ceil(self._zmin) 
if zmax is None: 

zmax = basics.custom.floor(self._zmax) 


for 


height in range(zmin, zmax + 1): 

for p in self.integer_points_in_bounding_rect_of_slice(height) 
if self.contains (p) : 
yield p 


def has_ld_greater_than(self, g) : 

"""Returns true iff ’self’ has lattice diameter greater than ’g’." 
for [p, q] in itertools.combinations(self.integer_points(), 2): 
if (p — q).gcd() > g: 
return True 
return False 


def _eq_(self, other): 

"""Returns true iff ’self’ and ’other’ represent the same 
polytope.""" 

return self.vertices == other.vertices 
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def _hash_(sel£) : 

return hash(self..vertices) 

def success!ve_minimum_at_most(self, value): 

"""Returns true iff the first successive minimum of ’self’ is at 
most ’value’.""" 

# Note that the first successive minimum of P ;= ’self’ is at most 

’value’ if and only if P — P contains a 

# non—zero point of 1/’value’ * Z''3; setting ’value’ =; a/b, this is 

equivalent to the fact that a * (P—P) 

# contains a non—zero integer point in b * Z''d 
value = Fraction(value) 

difference_body = PolytopeSd([(v — w) * value.numerator for v in 
self . .vertices for w in self..vertices 

if V != w]) 

for p in difference.body.integer.points(): 
if p.gcdO >= value . denominator : 
return True 
return False 

def vertices(self) : 

"""Returns a list of all vertices of ’self’.""" 
return list(self . .vertices) 


def _mul..(self, multiple): 

"""For a scalar ’multiple’, this returns ’multiple’ * ’self’.""' 
V=[v * multiple for v in self..vertices] 
return Polytopeld(V) 


def homo.copy(self, scaling, center): 

"""Returns 'scaling' * ’self’ + (1 — ’scaling’) * ’center’.""" 

V = [v * scaling + center * (1 — scaling) for v in self..vertices] 
return PolytopeSd(V) 


def all.facets.blocked(self): 

"""Returns true iff every facet of ’self’ contains an integer point 
in its relative interior.""" 
intpoints = list(self.integer.points()) 
for f in self..facets : 

if not self.is.facet.blocked(f, intpoints): 
return False 
return True 


def is.facet.blocked(self, facet, points): 
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"Returns true iff 'facet’ contains one point of 'points’ in its 
relative interior.""" 

227 for p in points: 

228 if facet.holds_with_equality(p) : 

229 if len([g for g in self..facets if 

g.holds.with.equality(p)]) == 1: 

230 return True 

231 return False 

232 

233 def is_l attice.free.extendable(self, points): 

234 """ 

235 Returns true iff there exists a point p in 'points’ \ 'self such 

that conv C self’ union {p}) is lattice—free. 

236 """ 

237 for p in points: 

238 if not self.contains(p) and Polytopeld(self.vertices() + 

[p]).is.lattice.free(): 

239 return True 

240 return False 

scripts.py 

1 """Contains very basic functions and objects needed for geometric 

computations. """ 

2 

3 import fractions 

4 

5 def gcd_three(a, b, c): 

6 """Returns the absolute value of the greatest common divisor of 'a', 

'b’, and 'c ’. """ 

7 return abs(fractions.gcd(a, fractions.gcd(b, c))) 

8 

9 def cu stom.floor(x): 

10 """Floor function that does not use floating point arithmetics when 

applied to fractions.""" 

11 return x // 1 

12 

13 def custom.ceil(x): 

14 """Ceiling function that does not use floating point arithmetics when 

applied to fractions.""" 

15 return — custom.floor (— x) 

16 

17 def normal_vector_of_triangle(points): 

18 """Returns a primitive nonzero vector 'v' in Z''3 that is perpendicular 

to the plane containing the 

19 (three-dimensional) points ’points[Q]’,’poins[1]’,’points[2] 
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20 


a = points [1] — points [0] 

21 b = points [2] — points [0] 

22 normal = a.cross(b) 

23 if normal.is.zero(): 

24 return Vector3d(0,0,0) 

25 return normal.primitive_vector() 

26 

27 def vol ume_of_simplex(vertices): 

28 """Computes the volume of the 3—dimensional simplex whose vertices are 

'vertices ’.""" 

29 assert len(vertices) == 4 

30 a = vertices[1] — vertices[0] 

31 b = vertices[2] — vertices[0] 

32 c = vertices[3] — vertices[0] 

33 return abs(fractions.Fraction(a.cross(b) * c, 6)) 

34 

35 # - 

36 

37 class Vector3d(object): 

38 """Straight—forward implementation of basic functions on vectors with 

three components.""" 

39 def __ init_(self, x, y, z): 

40 """Creates a vector with coordinates C'x', 'y', ’z 

41 self.x = X 

42 self.y = y 

43 self.z = z 

44 

45 def _str__(self): 

46 return "[" + str(self.x) + ", " + str(self.y) + ", " + str(self.z) + 

47 

48 def_ repr_(self): 

49 return str(self) 

50 

51 def __eq__(self, other): 

52 return self.x == other.x and self.y == other.y and self.z == other.z 

53 

54 def __ne__(self, other): 

55 return self.x != other.x or self.y != other.y or self.z != other.z 

56 

57 def _hash_(self): 

58 return hash(str(self)) 

59 

60 def add__(self, other): 

61 return Vector3d(self.x + other.x, self.y + other.y, self.z + other.z) 
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63 


64 

65 

66 
67 


def _sub__(sel£, other): 

return VectorSd(self.x — other.x, self.y — other.y, self.z — other.z) 
def _neg__(self): 

return Vector3d(—self.x, —self.y, —self.z) 
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71 


72 

73 

74 


75 

76 


78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 


def _mul__(sel£, other): 

'"'"If ’other' is a scalar, this returns the vector that arises from 
’self’ by multiplying each component with 
’other’. If ’other’ is a vector, this returns the scalar product of 
’self’ and ’other’."”” 
if type(other) is VectorSd: 

return self.x * other.x + self.y * other.y + self.z * other.z 
return VectorSd(self.x * other, self.y * other, self.z * other) 


def _div__(self, other): 

"""Returns the vector that arises from 
component by 'other'.""" 
return Vector3d(self.x / other, self.y 


'self by dividing each 
/ other, self.z / other) 


def gcd(self): 

"""Returns the greatest common divisor of the entries of 'self'. 
return gcd.three(self.x, self.y, self.z) 


def is_ zero(self): 

"""Returns true iff all entries of 'self are zero." 
return self.x == 0 and self.y == 0 and self.z == 0 


def is_ integral(self): 

"""Returns true iff all entries of 'self are integer numbers.""" 
return all([type(t) is int for t in [self.x, self.y, self.z]]) 

def primitive_vector(self): 

"""Returns an integral vector that is parallel to 'self and whose 
entries have gcd one.""" 

X = fractions.Fraction(self.x) 
y = fractions.Fraction(self.y) 
z = fractions.Fraction(self.z) 

V = VectorSd(x.numerator * y.denominator * z.denominator, \ 

y. numerator * x.denominator * z.denominator, \ 

z. numerator * x.denominator * y.denominator) 

return v / v.gcdO 
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103 

104 

105 

106 

107 

108 


def 


cross(self, other): 
"""Returns the 3—dimen 
return Vector3d(self.y 
self.z 
self.X 


sional cross product of 'self 

* other.z — self.z * other.y, 

* other.X— self.x * other.z, 

* other.y— self.y * other.x) 


and 


’other ’ 
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